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ABSTRACT 

We propose a new mechanism for the prompt emission of gamma-ray burst. 
In our model electrons are continuously accelerated in the post shock region via 
plasma turbulence. Using the Monte Carlo technique, we mimic the second-order 
Fermi acceleration due to plasma turbulence and obtain photon spectra. Since 
the acceleration balances with the synchrotron cooling, the observed low-energy 
spectral index is naturally explained. The resultant spectra can be consistent 
with observed spectra at least below ~ 1 MeV. The model also predicts delayed 
GeV-TeV emission due to inverse Compton and broad pulse profile of optical 
emission in some cases. Although nontrivial assumptions are required to repro- 
duce MeV-GeV power-law spectra, the model implies the possibility to explain 
various kinds of luminosity correlations. 

Subject headings: cosmic rays — gamma rays: bursts — gamma rays: theory - 
radiation mechanisms: 



Introduction 



In the wid ely discussed internal shock scenario (see, e.g., reviews by iPiran 1 12005 



Meszaros 1 120061 ). the prompt emission of gamma-ray bursts (GRBs) is due to collisions 
among inhomogeneities within ultrarelativistic outflows, which lead to formation of shocks. 
The nonthermal photons, whose typical energy e p ~ a few hundred keV, are emitted from 
shock-accelerated electrons in highly magnetized plasma. However, several open problems 
for the internal shock model have been pointed out such as the radiative efficiency, various 
kinds of luminosity correlations, and so on. In this paper, we focus on the two open prob- 
lems: the energy transfer problem and the low-energy spectral index problem. The standard 
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model postulates that a large fraction of the kinetic energy carried by protons should be effi- 
ciently converted into that of relativistic electrons. However, it is apparent that the Coulomb 
interaction cannot transport the internal energy of heated protons into electrons to achieve 
energy equipartition, because the timescale of the Coulomb interaction is much longer than 
the dynamical timescale. While the simple first-order Fermi acceleration at the shock front is 
assumed to transfer the energy into electrons in the standard model, some unknown plasma 
processes may play an important role in the energy transfer. 

The second problem is in the spectral shape of the prompt emis sion. The observe d 



spectra of GRBs are well fitted with the conventional Band function flBand et al. Ill993l ); 



the photon number spectrum oc e a exp [—(2 + a)e/e p ] below (a — f3)e v /{2 + a), and oc e 13 
abov e it. The typical fitted value of the low energy spectral index is a = —1.0 (IPreece et al. 
2000l ). while the standard model predicts that photons from cooled electrons dominate 



the low energy region below e p , which leads to a = —1.5. To resolve this problem sev- 



eral alternative models , such as the thermal emission from photosphere (IMeszaros fc Rees 



2 OOOl ; lloka et al. 1 120071 . and references therein), have been considered. The Klein- Nishina 
effect on synchrotron self-Compton (SSC) process, which can affect the l ow energy syn- 



chrot r on spectrum, has been discussed freq uently (IDerishev et al. 



20091 ; iNakar et al. 1 12009 : IWang et al. I |2009j). Recently, IPe'er fc Zhang 



2001 



Bosniak et al. 



J2006h have sug- 



gested that the decay of magnetic fields (§|2]) may resolve the problem in the low energy 
spectral index. However, in that model, we have not yet found a reason why the decay 
timescale should always be comparable to the cooling timescale In a particle-in-cell 

(PIC) simulation for elect ron-positron pla s ma (IChang et al. 1 120081 ) . the decay timescale is 
close to the requirement of IPe'er fc Zhang I (120061 ). but long-term evolution of magnetic fields 
is still contro versial, partially bec ause the effect of accelerated particles is not quantitatively 



unveiled yet (IKeshet et al. Il2009i ). 



Motivated by these problems, we propose an alternative model for the prompt emission 
of GRBs (SJ). In our model electron-heating (second-order Fermi acceleration) due to plasm a 



turbulence continues during photon emission in shocked plasmas (jGhisellini fc Celotti 



19991 ) 



so that the resultant spectral index of low energy photons can be consistent with the obser- 
vations at least below ~ 1 MeV. Since the assumed heating timescale is longer than that in 
the standard scenario, we call this model the "slow heating model". After the free energy 
for the plasma instabilities is dissipated, the magnetic fields may decay and the synchrotron 
emission will cease. The further possibility to reproduce MeV-GeV power-law spectra within 
the framework of this model is discussed in §3 In addition, our model naturally predicts 
delayed GeV-TeV emission due to inverse Compton (IC) and broad pulse profile of optical 
emission under certain conditions (§E}. Finally, we summarize the results of our model, 
and compare them with the observed luminosity correlations in GRBs (^7J). We should 
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note that our syn chrotron model is differe nt from the IC m odels using the same terminology 
"slow heating" in lPe'er et al. I (120061 ) and lGiannios I (120081 ). in wh ich the electrons are heated 



slowly in a timescale comparable to the shell-expans ion timescale (jGhisellini fc Celotti 111999 



Stern fc Poutanen 112004 : IVurm fc Poutanen 1120091 ). 



2. Magnetic Field: Generation and Decay 



In the standard scenario, the magnetic fields are assumed to be generated/amplified 
in the region around s hocks via plasma instabi l ities, such as the Weibel instability (e.g. 



Kato 



Kazimura et al. 1998: Medvedev fc Loeb 1999: Silva et al. 2003: Nishikawa et al. 2005 



2007h or the tw o-stream instability due to high-energy particles accelerated at the 



shock (e.g., iBell II2004I ). The difference in temperatures of electrons and protons may also 
arouse some plasma instabilities. Generated magnetic fields efficiently interact with particles, 
whose Larmor radii are comparable to the typical scale of the turbulence. This situation is 
definitely different from the ideal magnetohydrodynamic (MHD) approximation. For exam- 
ple, where the electron Larmor radii are finite, the off-diagonal t erms in the electron pressur e 



tensor could appear and catalyze magnetic reconnection (e.g.. iMozer and Pritchett II2009I ). 
The energy of magnetosonic pertur bations can be transferred to resonant particles via the 
transit-time damping process (e.g.. ISchlickeiser fc Miller Ill998l ). It is natural, therefore, to 



consider that the generated magnetic fields may decay via interaction with particles after 
the free energy for instability excitation (anisotropy, inhomogeneity, different temperatures 
of electrons and ions etc.) is dissipated with a decay timescale tdcc- Actually, recent PIC sim- 
ulations of electron-po sitron plasmas show that magnetic turbu lences induced by the Weibel 
instability decay (e.g.. lChang et al. 1120081 ; iKeshet et al. 1 120091 ) . For t > t dec , electrons stop 
emitting photons via synchrotron radiation. Therefore, the decay of magnetic fields may 
suppress the photon emission from cooled electrons, which resolves the problem in the index 
a. 



Numerical Model: Standard Case 



First let us revisit the effects of decay of magnetic fields with the standard manner of 
the electron injection. In Figure [U, changing the decay ti mescale, we show GRB spectra 
obtained by numerical calculations with the same code in lAsano fc Inoue I (120071 . details 
will be explained in the following section). Throughout this paper, all spectra are shown in 
terms of the observed fluence versus photon energy, assuming a GRB redshift of z = 0.1. 
The vertical axes denote ef(e), so that photon spectra with a spectral index a are plotted as 
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oc e a+2 . The model parameters are estimated as follows. The emitting region for a pulse is a 
homogeneous shell expanding with the Lorentz factor T at radius R from the central engine. 
We adopt I = R/T for the co moving width of the shell, so that the pulse timescale in the 



observer frame is At = R/Y 2 c (jSari fc Piran 1119971 ). Here, we choose parameters, V = 300 



At = 0.1 s, which implies R = 2.7 x 10 14 cm. The energy density of accelerated electrons in 
the shell U c = e c U (U is the total energy density of the shocked plasma) is a parameter that 
can be directly related to the isotropic-equivalent energy of photons from a single pulse E sh 
(here we adopt 10 51 erg) as E s ^ = U e V, where V = 4nR 3 /T is the comoving volume. 

In the standard scenario, relativistic electrons are injected at the shock front with a 
power-law energy distribution iV(7 e ) oc % p for 7 e > 7 e>m , where 7 e is the electron Lorentz 
factor in the plasma rest frame. This scenario requires a sharp low-energy cutoff for the 
electron injection spectrum; the minimum Lorentz factor 7 e m is evaluated in the literature 
by giving the energy density of electrons U e = e e U together with the total number density 
of electrons. Therefore, in the standard scenario, 7e, m has been conventionally described by 
the phenomenological parameter e e , though the energy scale corresponding to 7e, m should 
be derived from physics in relativistic plasmas. Here, instead of e e , we take 7 e , m to be a 
parameter, because we do not concern the non-observable parameter U. 

The photon energy e p corresponding to 7e )m is given by 

e p ~ :-r. i 

m e c 

The cooling timescale for electrons of 7 ejiri is written as 

, , v 6nm e c 

M7e,m) = ^ , (2) 

where <tt is the Thomson cross section. With a non-dimensional parameter e#, the magnetic 
energy density Ub = B 2 /8tt is given as e B U = (e B /e e )U e . In Figured! we set B = 3200 G 
and 7e jm = 3900, which correspond to es/e c = 0.1 and e p ~ a few hundred keV, respectively. 

We numerically follow electron cooling via synchrotron and IC emissions, adopting the 
Klein-Nishina cross section, and artificially stop the calculation after tdec to mimic the decay 
of magnetic fields. The effects of 77 pair production and synchrotron self-absorption are also 
taken into account. Since 7 e)m £ p /r > m e c 2 in our choice, the Klein-Nishina effect cannot be 
neglected for IC emission. The dynamical timescale tdyn = l/c = 30 s is much longer than 
the cooling time t c = 0.02 s. 



The blue curve with t dec = 0.01t c t c corresponds to the slow cooling case (jSari et al 



19981 ). In this case, e p is determined by the lowest energy of electrons that can cool within 
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Fig. 1, 

the decay timescale tdec- The assumed parameters are p = 2.5, E s h = 10 51 erg, T = 300 



Photon spectra for the standard model with decaying magnetic fields varying 

erg, 



R = 2.7x 10 14 cm (At = 0.1 s), B = 3200 G, and 7 e 



3900. 



the timescale tdec- As is well known, the index in the slow cooling case is a = — {p + l)/2, 
which is softer than the typical observed a for our choice of p — 2.5. If we adopt a very hard 
injection index p ~ 1, a can be ~ —1, but the high-energy index j3 = —(p + 2)/2 ~ —1.5 
contradicts the typical value f3 < —2. On the other hand, for the two cases, t& ec = 0.1t c 
and t c (green and blac k curves), the index becomes a ~ — 1 below e p . These cases are what 



Pe'er fc Zhang I (120061 ) suggested to solve the problem of the low-energy spectral index. 



It is further seen in Figure [T] that for the case of tdec = 10t c (red curve) the spectrum 
shows a ~ —1.5. This is the prediction by the standard model as referred in the introduction: 
electrons injected with 7 e = 7 e m are cooled after t = t c , and the low-energy spectrum becomes 
soft owing to emissions from such cooled electrons. It is noted that this case shows a spectral 
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bump in the GeV band due to IC emission, whose contribution is boosted up by enhancement 
of low-energy seed photons. We can also see cutoffs above 10 GeV and below 30 eV. They 
are 77-absorption and synchrotron self-absorption, respectively. 

The above results indicate that only with the case, tdec ~ £ c (7e,m), the decaying magnetic 
field can explain the low-energy index a. However, there is no definite physical reason to 



expect such a m atching between t dec and t c (7e,m)- The jitter radiation (IMedvedev 1 12000 



Fleishman 1 120061 ) instead of the synchrotron radiation is worthwhile to consider, because 
the typical scale of turbulence excited by plasma instabilities can be much shorter than the 
Larmor radii of radiating electrons. While the typical photon energy in the jitter radiation, 
which is determined by the coherence scale of the disturbed magnetic field, differs from in 
the usual synchrotron radiation, the introduction of the jitter radiation does not significantly 
change the low-energy spectral shape: a remains ~ —1.5 as long as tdec 3> (fast cooling). 



4. Slow Heating Model: start 

As we mentioned in §2], turbulent magnetic fields may be generated in the plasmas 
around shocks. Such turbulent waves may play a role in energy transfer from protons to 
electrons until the magnetic fields decay. In this section we present our new model, the 
slow heating model, to resolve the index problem. While the standard picture postulates a 
prompt acceleration of electrons, whose timescale is much shorter than £ c (7e,m), our model 
assumes slower energy transfer from the background plasma to electrons via some unknown 
plasma instabilities (see Figure [2]). 

In order to mimic the energy transfer we consider the second-order Fermi acceleration, 

eration 
et al. II2004T) or co- 



even though there may exist not only Alfven waves but also other types of acce 



mechanisms such as electric fields around ion current channels (IHededa 



herent wave-particle interactions resulting from parametric instabilities (IMatsukiyo fc Hada 



20091 ). which can also contribute to particle acceleration. When a particle is scattered by 
a wave or magnetized cloud preserving its energy in the wave (cloud) frame, twice Lorentz 
transformations give us energy gain due to this collision £ = AE/E as 

C = 7o (1 - Pom + /V2 - ~ 1, (3) 



where 70 = 1/yl — /3q, /ii, and fi' 2 are the Lorentz factor of the wave (cloud), cosines 
of incident angle in the reference frame, and scattering angle in the wave (cloud) frame, 
respectively. If the wave velocity is non-relativistic (@ < 1), the mean energy gain £ ~ (3% 
under the assumption of isotropic wave distribution JT[ = — f3 /2 and isotropic scattering 
/4 = 0. So the second-order Fermi acceleration is a slower acceleration process than the 
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Fig. 2. — Schematic pictures of the standard model and slow heating model. In the standard 
picture, radiation from the area behind the emission region cannot be neglected, as long as 
the magnetic field remains. On the other hand, the end of heating due to plasma turbulence 
and decay of the magnetic field are naturally simultaneous in the slow heating model. 



first-order one (£ ~ /3q) in non-relativistic cases. However, GRB internal shock is relativistic 
so that we can expect turbulent magnetic fields with /3q ~ 1. 



The Fokker-Planck equation for ultrarelativistic particles can be written as 

D EE 



ON 
~dt 



d_ ON 
dE EE dE 



d_ 

dE 



E 



E 



cool 



N 



(4) 



where Dee is the energy diffusion coefficient (see, e.g. JLiu et al. II2006I ). Defining the mean 
free time of particles t co u, we can write 



(5) 
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and the acceleration timescale t acc = £ C oii/£- At present we have no reliable model of relativis- 
tic turbulence in GRBs. For reference, let us look in stochastic acceleration in non-relativistic 
plasma. When we express the diffusion coefficient as D EE oc E n (t acc oc E 2 ~ n ), we obtain 
n = m for isotropic Alfven turbulence of sp ectral energy density per unit wavenumber 
Wk oc k~ m (see, e.g., iMiller fc Ramaty I Il989l). The model with n = 2 is often adopted 



for small scale MHD turbulences (see, e.g., iLiu et al. 1 120061 ) . Another value n = 5/3 is 



frequently used for the Kolmogorov turbulence. The strong turbulence limit (Bohm limit), 
where the mean free path becomes comparable to Larmor radii, corresponds to the case 
of n = 1, where the depend ence on m disappea rs. Relativistic shocks in electron-positron 
plasmas in PIC simulations (jChang et al. 1 120081 ) generate magnetic turbulence with m = 
for small k (kc <C plasma frequency), but m ~ 2 for large k. Although we have no definite 
shape of Dee for GRBs yet, future long-term PIC simulations may reveal the property and 
evolution of magnetic turbulence. 



From eq. (jl]) we may write AE 2 = (£ - 
function of £ per collision Gaussian form, 



?)E* 



so that we assume the probability 



p(0 



2na 



exp 



(g-0 s 

2a 2 



a 



(6) 



Hereafter, we assume a constant value of £ = 0.1. Considering synchrotron and IC emissions, 
the energy loss rate due to radiative cooling is expressed as 



E, 



4 



cool 



cttc^Ub 1 + 



P h 



U 



B 



(7) 



where U p h and K{j c ) are the photon energy density and the correction coefficient due to the 
Klein-Nishina effect, respectively. In the Thomson limit, K{p(^) = 1. 



We employ the Monte Carlo numerical code of lAsano fc Inoue I (120071 ) to follow the 
radiative cooling and stochastic energy gain/loss processes according to eqs. ([7]) and (jSJ) 
with a time step, 

6t = min(t coll /30, E/E cool /30, t dec /30), (8) 

and at t — tdec, we artificially halt the calculations to mimic the decay of magnetic fields. 
For each time step, we judge the occurrence of collision and estimate energy loss due to 
radiation using random numbers. If a collision occurs, the energy gain/loss due to the 
collision is counted with evaluated £. 

Since highl y disturbed magnetic fields are ass umed, it is meaningful to consider the 
jitter radiation ( Medvedev 2000 ; Fleishman"! 2006 ). But, for simplicity, we consider usual 



- 9- 



synchrotron radiation, using the synchrotron function (IRybicki and Lightman I Il979l ). As 
for IC emission we numerically estimate the spectral photon emission rate and i^(7 e ) by 
integrating t he photon energy distributio n given in advance with the Klein-Nishina cross 



section cx^n 

jRvbicki and Lightman Ill979l ). We assume a uniform and isotropic photon field 



within a shell with width I = R/T in the shell frame. To obtain photon spectra, the energy 
distributions of photons and particles are simulated iteratively until the resultant spectrum 
and presupposed spectrum are identical. 

In addition we take into account 77 pair production and synchrotron self-absorption. 
However, these processes are not so import ant in this paper , so th at we omit the explanation 



of the method to include these effects (see lAsano fc Inoue 1 120071 ). 



The studies for plasma turbulences in the post shock region by many authors are on- 
going now. Alth ough remarkable de velopment is seen in recent PIC simulations and MHD 



simulations (e.g. JZhang et al. 1120091 ). a definite picture of shocked plasma is not understood 
yet. Here, we consider a simple toy model assuming that Dee E 2 , which means that t co n 
does not depend on the energy of electrons. Although we take into account IC emission, 
it is not a main subject to discuss in this paper. In order to concentrate on synchrotron 
photon spectra, we adopt a stronger magnetic field B = 10 4 G. The other parameters are 
the same as those in £J3] except for the electron injection. The typical photon energy is 
expected to be emitted from electrons, whose energy loss rate is balanced with the second- 
order Fermi acceleration. Therefore, we adjust £ coU to make t acc = t c at 7 e = 7 typ = 3100 
(t c (7t yp ) = ^c,typ ~ 2 x 10~ 3 s) that implies the typical photon energy £ p (7 typ ) ~ a few 
hundred keV. The Klein-Nishina effect is important for electrons of 7 e = 7 typ even in this 
case. The number of electrons is roughly adjusted to make = 10 50 -10 51 erg considering 
the heating rate and tdec (we may not exactly forecast the final photon energy in advance). 
Below (above) 7 typ the acceleration timescale is shorter (longer) than the synchrotron cooling 
timescale. The heating due to turbulence reduces the effective number of electrons below 
7t yp so that the low-energy photon spectrum is expected to be harder than the standard 
one (the red line in Figured]) even for tdcc £c,typ- At t = electrons are injected with 
monochromatic energy of 7 e = 7 in j = 7 typ /10. We have confirmed that a run with 5000 par- 
ticle histories is enough to converge. In order to verify that the low-energy photon spectra 
become hard enough (a ~ —1) even for a longer decay timescale than the cooling timescale, 
we adopt tdec = 30t Cityp . The result is shown in Figure [3j 

Photons at the spectral peak are emitted from electrons of % ~ 7 typ as anticipated in 
advance. Since electrons are accelerated immediately, the injection parameter 7 in j does not 
affect the resultant spectrum very much. The low-energy spectral index is well approximated 
as ~ —1. The spectral bump at ~ GeV is due to IC emission, whose contribution is small 
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because of high B (es/e c ~ 3) and Klein-Nishina effect. Compared to Figured], the lower 
density of target photons (> 2 MeV) weakens the 77 absorption effect on IC emission. The 
overall shape of the spectrum is different from the Band function. We magnify the spectrum 
for 10 keV-1 MeV range in Figure E^b). In this energy range, the spectrum with artificial 
errors (10%) does not contradict the Band function very much, even though the model is 
quite simple. 



5. Slow Heating Model: modification 



While the model spectra in §3]may be fitted with the Band function below 
GRBs show power-law spectra in the MeV-GeV range with (3 ~ 



MeV, some 



sec 



Abdo et al. Il2009al . 



as one of the recent examples). As numerous simulations for the broadband prompt emis 



2007 



sion s pectrum have shown (e.g., iPe'er et al. 1120061 ; iGupta fc Zhang 1 120071 ; lAsano fc Inoue 



Bosnjak et al. 1120091 ). the usual simple power- law injection for electrons can easily 



reproduce the MeV-GeV power-law spectra. In the slow heating model, one of the simplest 
interpretation to overcome this difficulty is that such power-law spectra are superpositions 
of multiple components with different e p . This explanation requires fine adjust of the ampli- 
tude of multiple components. Although this interpretation remains viable so far, we search 
for alternative ideas within the slow heating model in this section. 

In §HJ the index n of energy dependence of _Dee(cx E n ) was taken 2. To make the 
resultant energy spectrum harder, we first test the cases with n > 2 (note that t acc < t c 
for n > 3 in higher energy range). Our simulations show that the spectral shape becomes 
close to a power-law function as n increases. However, even for an extreme choice of n = 3 
(t acc oc E^ 1 ), the resultant index (3 = —3.2 is still not hard enough. 

Since we consider the decay of the magnetic fields, it is natural to let Dee = 2£E 2 /2t co n 
also time (or equivalently distance from the shock front) dependent. Our toy model assumes 
a power-law shape as t co \\ oc 7°t x with upper and lower limits. Namely, 



t. 



coll 



mm 



c,typj 



max < t 



j ^min 



(9) 



where t m \ n = £t Ci t yp /100 in our model. In this case, the initial short timescale of acceleration 
makes e p higher, and as the acceleration timescale elongates with time, e p will be settled 
around a few hundred keV. Here, to harden spectra, we adopt 7 in j = 107 typ . Since the 
acceleration time from 7 in j to 1007 typ is ~ 10t m i n ^> t m i n , the highest photon energy may 
be ~ 100 MeV emitted from electrons of 7 e = 7 m j (10 GeV photons from electrons of 
7 e = 1007t yp may not be produced so much). Assuming t dec = 300t Cityp (other parameters 
are the same as before), we calculate spectra (see Figure @J. 
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We can see that the spectra above e p for x — 0.3 and 0.4 are well approximated by 
power-law functions. Even for x — 0.5, the spectrum from ~ 500 keV to ~ 3 MeV can be 
accepted as a power-law function. The peak energy e p for x — 0.3 becomes above MeV, 
because t acc at t = t^ ec is still shorter than the final acceleration timescale assumed in 



advance, ££ c ,ty P - For x — 0.5, after t = 10t Ci t yp <C td, 



300t, 



c,typ, 



the timescale t co ii attains 



the upper limit ££ c ,ty P , which makes 7 e ~ 7t yp . Compared to the timescale of stay around 
7; n j, the longer stayover around 7 typ yields the spectrum bump around a few hundred keV 
as set in advance. So one may easily understand that the spectra in this model depend on 
the timescale tdec- Figure [5] shows that the spectral bump around e p grows as tdec extends. 
It is apparent that the upper limit for t co n, whom we set up to adjust e p , causes the bumps. 
If t co ii is elongated monotonically, e p is determined by t co ii at t = tdec as the case of x — 0.3 
in Figure HI 

The models of Dee in this section are toy models to demonstrate the capability of 
MeV-GeV power-law spectrum in the slow heating model. In addition, we have assumed the 
monochromatic injection of electrons. The actual plasma turbulence and electron injection 
mechanism may be more complicated than the models we tested. The strength of magnetic 
fields may evolve with -Dee, or the first-order Fermi acceleration or surfing/ drift /wake-field 



acceleration (see, e.g., lAmano &: Hoshino 1 120071 ; iHoshino 1120081 ) may work as the injection 



mechanism at shock front. Although we need the nontrivial shape of Dee to reproduce 
the high-energy power-law spectra, the actual GRB plasmas may provide favorable condi- 
tions for Me V-GeV emissio n s. PI C simulations can be strong tools to verify this scenario. 
For example, IChang et al. I (120081 ) show that the energy density of magnetic turbulence in 
electron-positron plasma evolves as oc t~ 2 / 3 initially, then steepens to oc t~ l later. Such 
results encourage the model we discussed in this section. 



6. Delayed Emission 

In the simulations discussed in the previous sections, we have artificially halted the 
calculations at t = tdec with constant magnetic fields. However, the actual magnetic fields 
may not disappear suddenly, and we may expect residual magnetic fields at t > tdec- In this 
section, we consider emissions after the decay of magnetic fields. Here, we assume a simple 
exponential decay and residual magnetic fields as 



J B = max( J B e- t /* d -, J B min ), 



(10) 



- 12 - 



where Bq and B m [ n are constants. As for the acceleration timescale, to get rid of the heating 
effect smoothly, we assume the rapid evolution of t co u as 

£ coD = min (^ c ,t yP e {t/ * dcc)2 , t dya ) . (11) 

We numerically follow electron cooling/heating and photon emission during a period of 
£sim = ^dyn = I / c tdoc with 7i n j = 7 typ /10 and tdec = 30t Cityp (see Figure [6]). Two parameter 
sets are adopted; one is the same as that in £JH (r = 300, 7 typ = 3100, R = 2.7 x 10 14 cm, 
B = 10 4 G), and another parameter set describes a higher T case with the same e p and At: 
r = 800, 7 typ = 7800, R = 1.9 x 10 15 cm, and B = 530 G. For the final magnetic fields, 
-Bmin = 1 G is adopted in both the two cases, though there is no clue to the residual magnetic 
fields at present. For t <C idee, the photon emission mechanism is the same as those in §HJ 
However, after t = idee, the main cooling mechanism is switched from synchrotron to IC, 
because the photon density is assumed to be constant within the shell of width /. 

Comparing the thin dotted line with solid lines in Figure El it is clearly shown that the 
residual energy of electrons at t — tdec is emitted via IC emission in GeV-TeV ranges. The 
typical photon energy due to IC largely depends on T. As shown by the long dashed line 
in Figure El 77 absorption affects the final spectrum for F = 300 above 10 GeV (close to 
the cases of Figured]), while it is negligible for T = 800 owing to the lower photon density. 
Therefore, this case indicates a delayed onset of GeV-TeV photons compared to MeV photons 
with a timescale of ~ t^ ec /F. On the other hand, the spectral shape in the low-energy region 
is not altered in this model as seen in Figure El 

Recent GRBs detected by Fermi-LAT tend to show su ch a delayed onset o f high-energy 
(> 100 MeV) emission (GRB 080916C, GRB 080825C, etc.: lAbdo etal~ll2009al Jbll However, 



it is noted that the broadband spectral shape for GRB 080916C is not well reproduced by 
our simple toy model. Future observations with Fermi or Cerenkov telescopes will testify 
the model prospects and yield a clue to the refinement of the model. 

Next, let us consider photon emission in longer timescales than the dynamical timescale 
td y n- In the standard scenario, emission after the dynamical timescale may be negligible 
because of the fast cooling of electrons. However, in our scenario, the residual electron 
energy may be released v ia synchrotron radiation, which can contrib ute to optical emissions 



as seen in some GRBs ( jVestrand et al. 1 120051 ; iBlake et al. 1 120051 ). Our results for such 



cases are shown in Figure El where the parameters are the same as those for F = 300 in 
Figure El but t sim = 10td yn , B min = 1, 10 and 30 G, and rapid decay of the photon density 
oc exp [— (t/tdec) 2 ]- For simplicity, we neglect the effects of adiabatic cooling due to shell 
expansion. 

As seen in Figure H, optical synchrotron emissions become more luminous than the flu- 
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ence of the power-law extrapolation from the X-ray spectra. Most of the optical photons are 
emitted after t = t djn in the weak magnetic fields so that broader pulse profile is expected 
for optical than 7/X-ray bands . Such longer variabi lity timescales in prompt optical emis- 



sion is seen in GRB 080319B (IRacusin et al. 1 120081 ) . though the flux is brighter than the 



extrapolation from the X-ray spectra by 3-4 orders of magnitude (1-2 orders in our results). 
More luminous optical emissions as compared to X-ray may be possible, if we change tdcc- 
Since our objective in this paper is not to reproduce spectra of specific GRBs, we do not 
further discuss the fraction of optical flux here. 



7. Summary and Discussion 

Motivated by the energy transfer and low-energy spectrum problems, we propose a new 
model to reproduce GRB prompt emissions. In this model, electrons are continuously heated 
via plasma turbulences within a timescale longer than the cooling timescale. The acceleration 
timescale is assumed to be much longer than that in the standard picture. Emissions from 
cooled electrons are suppressed so that the low-energy spectral index a is close to the observed 
value —1. At least below MeV, the model spectrum does not contradict the Band function 
very much. Considering that most of the GRB spectra were obtained in energy ranges below 
MeV, the model spectra may be consistent with a large fraction of GRBs. In order to explain 
power-law spectra in MeV-GeV range observed in some GRBs, we need a superposition of 
multiple components with different e p , or nontrivial shape and evolution of the diffusion 
coefficient. Our model, under certain conditions, predicts delayed GeV-TeV emission via 
IC or delayed optical emission with broad pulse profile via synchrotron. We expect that 
the accumulation of many GRB observations will verify the characteristics predicted by our 
model in near future. 

Roughly speaking, the energy release in this model is estimated as T N(t dec /t c ,ty P )ltypm e c 2 , 
where N is the total number of accelerated electrons. Given the total isotropic energy 
E iso = 4nIi 3 U, the required number of electrons depends on t dec . Thus, the fraction of 
accelerated electrons can be much less than unity in this model, while many autho rs have 



frequently assumed that all electrons are accelerated (see lEichler fc Waxman II2005I ). 



One interesting point in our model is that the physical explanation for the spectral peak 
energy e p is clear. The balance between synchrotron cooling and heating provides us the 
typical electron energy 7t yp m e c 2 , from which we can estimate e p . The spectral peak and low- 
energy index are reproduced by this mechanism unless t acc 3> t c . One remaining problem is 
why shocked plasma in GRBs always adjust t co ii to make e p observed range of ~ 1 MeV. To 
explain this we need another assumption; for example, D EE oc Ub, which implies t acc oc B~ 2 . 
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The balance between cooling time t c oc B 2 and t acc gives us 7 typ oc B°. Since we may write 
the luminosity as L iso = E iso / (R/cT 2 ) = 47r(e c / € B )UbR 2 T 2 , we obtain 

e p oc r£ 7t 2 yp oc (e c /e B )-V 2 LH 2 R-\ (12) 

where Inde pendence disappears. T hese results are consistent with the Yonetoku relation 
(e p oc L-j^; lYonetoku et al. 1120041 ). except for the factor Rr 1 . Of course, the above as- 



sumption is not trivial. When we assume Dee oc [/§, the correlation becomes e p oc 
T 2x ~ 2 R 2x ~ 3 . If the jitter radiation is applicable, in which the typical photon en- 
ergy depends on the coherent scale of the field, the Yonetoku relation implies some corre- 
lations between the typical scale of turbulence and luminosity. In any case, fundamental 
studies of the long-term evolution of relativistic plasmas, based on PIC simulations, etc., are 
indispensable to verify the model and luminosity correlations. 



For the Fermi-L AT GRB 0809 16C, by a process of elimination, IZhang fc Pe'erl (120091 ) 
conclude that this GRB is emitted from a magnetically dominated outflow. In such a case, 
the dissipation of the bulk kinetic energy may not be due to internal shocks, because the 
Alfven velocity is very close to the light speed. In order to produce non-thermal particles, 
the dissipation of magnetic fields such as magnetic reconnection, etc. should occur in the 
outflow. The dissipation processes of magnetic fields may generate magnetic turbulences so 
that we can expect the second-order Fermi acceleration in this case too. 



We appreciate the anonymous referee for the useful advice. This work is partially 
supported by the Grant-in-Aid for Scientific Research, No. 21540259 from the MEXT of 
Japan. 
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Fig. 3. — (a) Spectrum for the slow heating model with constant £ and t co \\. (b) Zoom up 
of the upper figure (the area enclosed by the rectangle) for 10 keV - 1 MeV range. The data 
points are calculated spectra with 10 % error putted in by hand, and the solid line is the 
Band function with e p = 200 keV, a = —1.15, and (5 = —2.3. 
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Fig. 4. — Spectra for the model of t co n oc t x [Dee oc t~ x , x = 0.3, 0.4, and 0.5), 7 in j = 107 typ , 
and t dcc = 300t City p. 



- 19 - 



ef(e) [erg/cm ] 




e[eV] 

Fig. 5. — Spectra for the model of D EE oc t~ 0A , 7 in j = 107t yp , and different t dcc (100, 300, 
and 1000 t c , typ ). 
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Fig. 6. — Spectra for the model of B oc e _t// * dcc , 7i n j = 7 typ /10, tdec = 30t Ci t yp and t sim = td yn , 
and different T (solid lines, see the text). The thin dotted line is the same as the spectrum 
in Figure [3j where t sim = t dec t dyn and B is constant. The long dashed line is the IC 
component for V = 300 without 77 absorption effect. 
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